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SUMMARY 

A  comprehensive  analytical  procedure  was  developed  for  predicting  faults  in  gear  transmission 
systems  under  normal  operating  conditions.  A  gear  tooth  fault  model  is  developed  to  simulate  the  effects 
of  pitting  and  wear  on  the  vibration  signal  under  normal  operating  conditions.  The  model  uses  changes  in 
the  gear  mesh  stiffness  to  simulate  the  effects  of  gear  tooth  faults.  The  overall  dynamics  of  the  gear  trans¬ 
mission  system  is  evaluated  by  coupling  the  dynamics  of  each  individual  gear-rotor  system  through  gear 
mesh  forces  generated  between  each  gear-rotor  system  and  the  bearing  forces  generated  between  the  rotor 
and  the  gearbox  structure.  The  predicted  results  were  compared  with  experimental  results  obtained  from  a 
spiral  bevel  gear  fatigue  test  rig  at  NASA  Lewis  Research  Center.  The  Wigner-Ville  Distribution  (WVD) 
was  used  to  give  a  comprehensive  comparison  of  the  predicted  and  experimental  results.  The  WVD 
method  applied  to  the  experimental  results  were  also  compared  to  other  fault  detection  techniques  to 
verify  the  WVD’s  ability  to  detect  the  pitting  damage,  and  to  determine  its  relative  performance.  Overall 
results  show  good  correlation  between  the  experimental  vibration  data  of  the  damaged  test  gear  and  the 
predicted  vibration  from  the  model  with  simulated  gear  tooth  pitting  damage.  Results  also  verified  that  the 
WVD  method  can  successfully  detect  and  locate  gear  tooth  wear  and  pitting  damage. 
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I.  INTRODUCTION 

In  the  last  two  decades,  the  use  of  gear  transmissions  in  both  defense  and  commercial  applications 
has  substantially  increased.  With  the  demand  for  higher  power  and  performance,  premature  failures  in 
transmissions  often  result  in  financial  losses,  and  sometimes  even  lead  to  catastrophic  consequences.  In 
the  aerospace  industry,  one  of  the  major  concerns  is  with  gear  fatigue  failures  in  rotorcraft  transmission 
systems.  Large  vibrations  in  gear  transmission  systems  usually  result  in  excessive  gear  tooth  wear  and 
possible  tooth  crack  formation  which,  in  turn,  leads  to  premature  gear  failure.  Thus,  it  is  important  to  un¬ 
derstand  the  dynamics  of  a  transmission  system  over  a  variety  of  fault  conditions,  as  well  as  under  nomi¬ 
nal  conditions.  With  this,  methods  can  be  explored  to  detect  and  assess  the  magnitude  of  the  gear  damage 


present.  Due  to  limitations  in  the  number  and  types  of  experiments  that  can  be  performed,  the  only  practi¬ 
cal  means  of  obtaining  this  type  of  data  is  through  analytical  simulations. 

The  major  objective  of  the  research  reported  herein  is  to  develop  and  verify  a  model  to  predict  the 
vibration  of  a  transmission  system  with  the  effects  of  gear  surface  pitting  and  wear.  To  simulate  the  vibra¬ 
tion  of  the  transmission  system,  the  equations  of  motion  were  established  individually  for  each  rotor-gear- 
bearing  system.  The  effects  of  tooth  wear  or  surface  pitting  are  simulated  by  changes  in  the  magnitude 
and  phase  of  the  mesh  stiffness.  These  localized  changes  in  the  gear  mesh  are  incorporated  into  each  gear- 
rotor  model  for  dynamic  simulation  [Choy  1991,1993  and  Kahraman  1990].  The  dynamics  of  each  gear- 
rotor  system  are  coupled  with  each  other  through  the  gear  mesh  interacting  forces  and  the  bearing  support 
forces.  The  global  vibrations  of  the  system  are  evaluated  by  solving  thetransient  dynamics  of  each  rotor 
system  simultaneously  with  the  vibration  of  the  casing.  In  order  to  minimize  the  computational  effort,  the 
number  of  degrees  of  freedom  of  the  system  are  reduced  by  using  a  modal  synthesis  procedure  [Choy 
1991,1993]. 

To  verify  the  analytical  model  with  experimental  data,  the  dynamics  of  a  single  spiral  bevel  pin¬ 
ion  with  various  degrees  of  gear  tooth  damage  was  simulated.  Results  from  the  model  were  compared  to 
experimental  results  using  a  joint  time-frequency  analysis  method.  This  approach  was  chosen  because  of 
the  large  amount  of  information  represented  in  the  joint  time-frequency  results  which  cannot  be  repre¬ 
sented  separately  in  either  the  time  domain  or  the  frequency  domain.  The  joint  time-frequency  analysis 
will  provide  an  instantaneous  frequency  spectmm  of  the  system  at  every  instant  of  the  revolution  of  the 
pinion  while  a  Fourier  Transform  can  only  provide  the  average  vibration  spectmm  of  the  signal  obtained 
during  one  complete  revolution.  In  other  words,  the  time-changing  spectral  density  from  the  joint  time- 
frequency  spectra  will  provide  information  concerning  the  frequency  distribution  concentrated  at  that  in¬ 
stant  around  the  excited  instantaneous  frequency  which  cannot  be  obtained  in  a  regular  vibration 
frequency  spectmm.  The  joint  time-frequency  analysis  approach  applies  the  Wigner-Ville  Distribution 
(WVD)  [Boashash  1987,  Claasen  1980,  and  Shin  1993]  on  the  time  vibration  signal  of  the  system.  Some 
success  has  been  achieved  in  applying  the  WVD  to  gear  transmission  systems  [Forrester  1990  and 
McFadden  1991]  to  recognize  faults  at  various  locations  of  the  gear.  Other  fault  detection  techniques, 
including  frequency  domain  analysis  [Randall  1982  and  Taylor  1980]  and  several  time  discriminant  meth¬ 
ods  such  as:  FMA  [Stewart  1977],  NAA*,  and  NBA*  [Zakrajsek  1993,1994  and  Decker  1994]  are  also  used 
to  compare  with  and  verify  the  WVD  approach. 

Based  on  results  of  this  study,  some  conclusions  are  made  on  the  ability  of  the  developed  model 
to  simulate  gear  tooth  surface  damage,  and  the  ability  of  the  WVD  method  to  detect  damage  sufficient  to 
verify  the  model. 


II.  ANALYTICAL  PROCEDURE 

The  dynamics  of  the  ith  individual  gear-shaft  system  can  be  evaluated  through  the  equations  of 
motion  for  the  vibrations  of  a  individual  rotor-bearing-gear  system  as  shown  in  figure  1  [Choy 
1991,1993],  given  in  matrix  form,  as 

[^]{^}  +  [Kspi]  =  {4(0}  +  {4(0}  {4(0}  (1 

where  [AT]  and  [K^]  are  respectively  the  mass  and  shaft  stiffness  matrices  of  the  rotor,  {W^}  is  the  general 
displacement  vector  of  the  ith  rotor  in  its  local  coordinate  system,  and,  {F^j(r)},  {F^f^f)},  and  {4f(0}  are 
respectively  the  force  vectors  acting  on  the  ith  rotor  system  due  to  bearing  forces,  gear  mesh  interactions, 
and  mass-imbalances.  In  this  model,  the  dynamics  between  the  gearbox  and  the  rotor  are  coupled  through 
the  bearing  forces,  which  are  evaluated  by  the  relative  motion  between  the  rotor  and  the  gearbox.  The 
interactions  between  each  individual  rotor  are  coupled  through  the  gear  forces  generated  by  the  relative 
motion  of  the  two  mating  gears  at  the  mesh  point. 
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(2) 


The  equations  of  motion  of  the  gearbox  with  p  rotor  systems  can  be  expressed  as 

i=l 

where  represents  the  coordinate  transformation  between  the  ith  rotor  and  the  gearbox. 

The  bearing  forces  for  the  ith  rotor  can  be  evaluated  as 

I'i.w} = [c«]({>i;}  -  [4l{»'d}) + KK{'^}  -  [tcllB;,}) 

where  [C^,]  and  are  respectively  the  damping  and  stiffness  of  the  bearing,  is  the  coordinate 
transformation  matrices  for  the  gearbox  with  respect  to  the  ith  rotor,  and  W^-  is  the  casing  displacements 
at  the  rotor  locations. 

The  gear  forces  generated  from  the  gear  mesh  interaction  [Choy  1991,1993]  can  be  written  as 

few}  •  feiWl  +  feiW} 

where  { F„(0 }  is  the  vector  containing  the  gear  forces  and  moments  resulting  from  the  relative  rotation 
between  the  two  mating  gears  and  is  the  vector  containing  gear  forces  and  moments  due  to  the 

translational  motion  between  the  two  gears  [Choy  1988  and  Boyd  1989]. 

In  order  to  calculate  the  transient/steady  state  dynamics  of  the  system,  all  the  coupled  rotor  and 
casing  equations  of  motion  have  to  be  solved  simultaneously.  To  minimize  the  computational  effort,  the 
modal  transformation  [Choy  1991,1993]  procedure  was  applied  to  reduce  the  degrees  of  freedom  of  the 
global  equations  of  motion.  Using  the  undamped  mode  shapes  of  the  rotor  system  and  the  undamped 
mode  shapes  of  the  gearbox,  the  rotor  and  the  gearbox  displacements  were  transformed  into  the  modal 
coordinates  to  reduce  the  number  of  degrees  of  freedom  in  the  system  to  minimize  the  computational  ef¬ 
fort.  The  modal  equations  of  motion  for  the  rotor  systems  and  the  gearbox  vibration  were  solved  simulta¬ 
neously  for  the  modal  accelerations.  A  numerical  integration  scheme  was  used  to  integrate  the 
accelerations  to  obtain  velocities  and  displacements  at  each  time  step  for  transient  calculations. 


m.  EXPERIMENTAL  STUDY 

Using  the  spiral  bevel  gear  fatigue  test  rig,  illustrated  in  figure  2,  the  resulting  fatigue  damage  on 
the  pinion  is  illustrated  in  figures  3  and  4.  The  primary  purpose  of  this  rig  is  to  study  the  effects  of  gear 
tooth  design,  gear  materials,  and  lubrication  types  on  the  fatigue  strength  of  aircraft  quality  gears 
[Handschuh  1992].  Because  spiral  bevel  gears  are  used  extensively  in  helicopter  transmissions  to  transfer 
power  between  nonparallel  intersecting  shafts,  the  use  of  this  fatigue  rig  for  diagnostic  studies  is  practical. 
Vibration  data  from  an  accelerometer  mounted  on  the  pinion  shaft  bearing  housing  was  captured  using  a 
personal  computer  with  an  analog  to  digital  conversion  board  and  antialiasing  filter.  The  12-tooth  test 
pinion,  and  the  36-tooth  gear  have:  0.5 141 -in.  diametral  pitch,  35°  spiral  angle,  1-in.  face  width,  90°  shaft 
angle,  and  22.5°  pressure  angle.  The  pinion  transmits  720  hp  at  nominal  speed  of  14  400  rpm.  The  test  rig 
was  stopped  several  times  for  gear  damage  inspection.  The  test  was  ended  at  17.79  operational  hours 
when  a  broken  portion  of  a  tooth  was  found  during  one  of  the  shutdowns. 


IV.  VIBRATION  SIGNATURE  ANALYSIS  PROCEDURE 

Three  major  methodologies:  (A)  the  joint  time-frequency  approach,  (B)  the  frequency  domain 
approach,  and  (C)  the  time  domain  techniques,  were  used  in  this  study.  The  following  is  a  description  of 
the  three  methodologies: 
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(A)  Joint  Time-Frequency  Technique 


To  examine  the  vibration  signal  in  a  joint  time-frequency  domain,  the  Wigner-Ville  method 
[Boashash  1987  and  Claasen  1980]  was  used  in  this  study.  While  the  Fast  Fourier  Transform  (FFT)  tech¬ 
nique  can  provide  the  spectral  contents  of  the  time  signal,  it  cannot  distinguish  time  phase  change  during 
a  complete  cycle  of  operation.  In  other  words,  it  assumes  that  the  time  signals  are  repeatable  for  each  time 
data  acquisition  window  without  considering  the  effects  of  any  magnitude  and  phase  changes  during  the 
sampling  period.  The  Wigner-Ville  distribution  will  provide  an  interactive  relationship  between  time  and 
frequency  during  the  period  of  the  time  data  window.  The  comprehensive  representation  of  the  vibration 
signal  using  the  WVD  method  is  the  primary  reason  that  it  was  used  to  compare  the  predicted  and  experi¬ 
mental  vibration  results.  The  WVD  (Wigner-Ville  Distribution),  in  a  discrete  form,  can  be  written  as: 

L 

W^{nT,f)  =  2T^x{nT +  iT)x*[nT  -  iT) ■  (5) 

i=-l 


where 

W(tJ)  the  Wigner-Ville  distribution  in  both  the  time  domain  t  and  frequency  domain/ 
x(t)  the  time  signal 

T  the  sampling  interval 

L  the  length  of  time  data  used  in  the  transform 

To  allow  sampling  at  the  Nyquist  rate  and  eliminate  the  concentration  of  energy  around  the  frequency 
origin  due  to  the  cross  product  between  negative  and  positive  frequency  [Boashash  1987  and  Claasen 
1980],  the  analytic  signal  was  used  in  evaluating  the  WVD.  The  analytic  signal  s(t)  is  defined  as 

s{t)  =  x(r)  +  jH[x{t)\  (6) 

where  H[xit)]  is  the  Hilbert  transform  of  x(t).  However,  an  alternative  approach  can  be  used  to  calculate 
the  analytic  signal  using  the  frequency  domain  definition.  The  analytic  signal  s(t)  can  be  evaluated  by 
calculating  the  FFT  of  the  time  signal  x(t),  then  setting  the  negative  frequency  spectrum  to  zero.  The  ana¬ 
lytic  signal  can  be  obtained  by  evaluating  the  inverse  FFT  of  the  spectrum. 

To  simplify  the  computational  effort,  the  WVD  can  be  evaluated  using  a  standard  FFT  algorithm.  Adopt¬ 
ing  the  convention  that  the  sampling  period  is  normalized  to  unity,  it  is  necessary  only  to  evaluate  the 
WVD  at  time  zero.  Hence 


W^{0J)  =  2 


(7) 


where  k(i)  =  s(i)s*(i).  Equation  (13)  can  be  evaluated  using  the  discrete  FFT  algorithm. 

In  order  to  avoid  a  repetition  in  the  time  domain  WVD,  a  weighting  function  [McFadden  1991] 
was  added  to  the  time  data  before  the  evaluation  process.  Such  a  process  may  decrease  the  resolution  of 
the  distribution,  but  it  will  eliminate  the  repetition  of  peaks  in  the  time  domain  and  the  interpretation  of 
the  result  is  substantially  easier. 
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(B)  Frequency  Domain  Technique 


The  frequency  spectrum  is  found  by  applying  a  Discrete  Fourier  Transform  (DFT)  on  the  time 
averaged  signal  x{t\  such  that  the  spectral  components  are 


;=0  ^  ^ 


-jlnik'] 


(8) 


where 

x{f)  time  domain  signal 

X(k)  frequency  domain  signal 

T  sampling  time  interval 

N  number  of  data  points 

The  frequency  components  were  examined  in  the  frequency  domain  and  compared  with  those  obtained  at 
various  stages  of  the  fault  development  in  the  spiral  bevel  pinion. 

(C)  Time  Domain  Techniques 

Four  different  time  domain  techniques  for  early  detection  of  gear  tooth  damage  were  used  in  this 
study.  All  of  the  time  domain  techniques  were  applied  to  the  vibration  signal  after  it  was  time  synchro¬ 
nously  averaged.  These  techniques  are:  FMA,  NAA  ,  and  NBA  .  They  are  defined  as  follows. 

FMA  [Stewart  1977]  is  an  isolated  fault  detection  parameter  for  the  gear  tooth,  and  is  given  by  the 
normalized  kurtosis,  of  the  resulting  difference  signal  as 


FMA  = 


N±(d,-df 

1=1 _ _ 

■  JV  f 

1(4 -3f 

-1=1 


(9) 


where 

d{t)  A{f)  -  R{t) 

mean  value  of  d{t) 

A{t)  original  time  synchronous  signal 

R{t)  regular  meshing  components  plus  their  first  order  side  bands 
N  total  number  of  data  points  in  the  time  signal 

NAA*  [Zakrajsek  1993,1994  and  Decker  1994]  is  a  general  gear  fault  detection  parameter,  with 
trending  capabilities.  A  residual  signal  is  constructed  by  removing  regular  meshing  components  from  the 
time  averaged  signal.  The  first  order  sidebands  stay  in  the  residual  signal  and  the  fourth  statistical  mo¬ 
ment  of  the  residual  signal  is  then  divided  by  the  averaged  variance  of  the  residual  signal,  raised  to  the 
second  power.  The  average  variance  is  the  mean  value  of  the  variance  of  all  previous  records  in  the  run 
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ensemble.  This  allows  M4  to  compare  the  current  gear  vibration  with  the  baseline  of  the  system  under 
nominal  conditions.  NAA  is  given  by  the  quasi-normalized  kurtosis  equation  shown  below: 


NAA*{M) 


(10) 


where 

r  residual  signal 

r-  mean  value  of  residual  signal 

N  total  number  of  data  points  in  one  time  record 

j  data  point  number  in  time  record 

j  time  record  number 

M  current  time  record  number  in  run  ensemble 

An  enhancement  to  this  parameter  is  given  by  NAA*,  in  which  the  value  of  the  averaged  variance  is 
“locked”  when  the  instantaneous  variance  exceeds  a  predetermined  value  [Decker  1994].  This  provides 
NAA  with  enhanced  trending  capabilities,  in  which  the  kurtosis  of  the  current  signal  is  compared  to  the 
variance  of  the  baseline  signal  under  nominal  conditions. 

NBA  is  another  parameter  similar  to  NAA  that  also  uses  the  quasi-normalized  kurtosis  given  in 
equation  (4).  The  major  difference  is  that  instead  of  using  a  residual  signal,  NBA  uses  the  envelope  of  the 
signal  bandpassed  about  the  mesh  frequency.  Again,  as  with  NAA  ,  NBA  is  an  enhancement  to  the  NBA 
parameter,  in  which  the  value  of  the  average  variance  is  “locked”  when  the  instantaneous  variance  ex¬ 
ceeds  a  predetermined  value.  The  equation  for  NBA  is  given  below: 


NBA*{M) 


Nj^(s,-7f 


1=1 


(11) 


and 

s{t)  =  magnitude  of  |fc(f)  +  i{//[l>(f)]}|  (12) 

where 

l){t)  time  averaged  signal  bandpassed  filtered  about  the  meshing  frequency 

H[b{ty]  the  Hilbert  Transform  of  b(t) 

N  total  number  of  data  points  in  one  time  record 

i  data  point  number  in  time  record 

j  time  record  number 

M  current  time  record  number  in  run  ensemble 
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V.  DISCUSSIONS  OF  RESULTS 


A  series  of  photos  showing  the  deterioration  of  the  pinion  teeth  at  various  stages  of  the  test  are 
shown  in  figures  3  and  4.  Figure  3(a)  shows  the  initiation  of  a  small  pit  on  one  of  the  pinion  teeth  during 
the  first  shutdown,  at  about  5.5  hr  into  the  test.  As  the  test  progressed,  the  rig  was  shut  down  seven  more 
times  to  examine  the  severity  of  the  pitting  and  its  relationship  with  the  corresponding  vibrations.  Figures 
3(b),  (c),  and  (d)  show  the  increase  of  the  damaged  area  at  the  pinion  tooth  as  the  elapsed  time  increased 
to  6.55,  8.55,  and  10.03  hr,  respectively.  Note  that  in  figure  4(a),  at  12  hr,  the  damage  of  the  pinion  tooth 
increased  to  75  percent  of  the  tooth  surface.  At  this  stage,  pitting  also  initiated  on  the  adjacent  tooth  and 
continued  to  grow  as  the  time  increased  to  14.53  hr  (fig.  4(b)).  At  16.16  hr,  the  damage  has  grown  to  three 
adjacent  teeth  as  shown  in  figure  4(c).  The  test  was  terminated  when  a  tooth  fracture  was  found  on  one  of 
the  three  heavily  pitted  teeth  at  17.79  hr,  as  shown  in  figure  4(d). 

Figures  5  and  6  show  the  WVD,  the  time  signal  (to  the  left  of  the  WVD),  and  frequency  spectra 
(below  the  WVD),  of  the  spiral  bevel  pinion  vibration  at  various  stages  of  damage,  corresponding  to  the 
photographs  given  in  figures  3  and  4.  Figure  5(a)  shows  the  vibration  signature  due  to  initial  pitting,  as 
shown  in  figure  3(a),  at  around  200°  from  the  triggering  point  of  the  gear  at  5.5  hr  of  running  time.  As 
time  progressed  to  6.55  hr,  more  advanced  pitting  is  occurring  as  shown  in  figure  3(b),  the  WVD  ampli¬ 
tude  increases  following  by  a  short  term  amplitude  and  phase  change  at  260°  as  illustrated  in  figure  5(b). 
This  phenomena  is  also  evident  in  the  frequency  spectrum  with  the  existence  of  sizable  sideband  compo¬ 
nents.  Due  to  the  speed  increase  (after  the  shutdown)  at  8.55  hr,  the  overall  WVD  amplitude  increases 
substantially  as  shown  in  figure  5(c).  At  the  ranning  time  of  10.03  hr,  as  seen  in  figure  3(d),  pitting  be¬ 
come  more  advanced  as  can  be  detected  by  the  initiation  of  a  cross  pattern  in  the  WVD  (fig.  5(d)).  The 
pitting  of  the  tooth  is  more  pronounced  at  12.03  hr,  as  shown  in  figure  4(a),  and  can  be  verified  by  the 
distinct  cross  pattern  in  the  WVD  shown  in  figure  6(a).  The  high  concentration  in  the  WVD  energy  at 
14.53  hr  combined  with  the  further  increase  in  sideband  amplitudes  in  figure  6(b)  shows  the  initiation  of 
the  pitting  process  on  the  neighboring  teeth,  as  illustrated  in  figure  4(b).  At  16.16  hr,  as  seen  in  figure 
6(c),  the  WVD  pattern  changes,  as  two  cross  patterns  start  to  appear,  showing  more  advanced  damage  to 
the  adjacent  teeth  as  illustrated  in  figure  4(c).  The  discontinuity  of  the  mesh  frequency  component  with 
time,  as  shown  by  the  WVD  in  figure  6(d),  is  mainly  due  to  the  short  term  phase  change  caused  by  the 
fractured  tooth  and  the  pitting  of  the  two  adjacent  teeth,  as  shown  in  figure  4(d).  This  phenomenon  is  also 
confirmed  by  the  substantial  increase  in  the  sideband  amplitude  in  the  frequency  spectrum  in  figure  6(d). 

In  order  to  simulate  the  dynamics  of  the  gear  system  with  gear  tooth  damage,  two  mesh  stiffness 
models,  with  combinations  of  stiffness  reduction  and  phase  delay,  were  introduced.  A  simple  stiffness 
reduction  and  phase  change  model,  shown  in  figure  7(a),  was  used  to  simulate  the  single  tooth  surface 
pitting  damage  at  12  hr,  as  shown  in  figure  4(a).  The  resulting  experimental  vibration  signature  due  to  this 
damage  is  given  in  figure  6(a).  The  mesh  stiffness  model  given  in  figure  7(b)  consisted  of  a  combination 
of  stiffness  and  phase  changes  over  three  consecutive  teeth,  to  approximate  the  multiple  tooth  pitting 
damage  at  17.79  hr,  as  seen  in  figure  4(d).  The  resulting  experimental  vibration  signature  due  to  this  ex¬ 
tensive  damage  is  given  in  figure  6(d).  Additional  frictional  effects  were  also  added  to  the  model  to  simu¬ 
late  the  roughness  of  the  tooth  surface  due  to  pitting.  The  simulated  vibration  signamres  of  the  pinion  gear 
are  given  in  figure  8.  Comparing  figures  8(a)  and  6(a),  for  the  single  tooth  damage  case  at  12  hr,  one  may 
notice  the  similarities  between  both  the  frequency  spectra  and  the  WVD  display.  The  unevenness  in  the 
experimental  time  signal  is  mainly  due  to  the  modulation  of  frequencies  due  to  other  excitations  in  the  test 
]*jg^  \^ttich  are  not  included  in  the  model.  For  the  tooth  break-off  case  at  17.79  hr,  both  the  numerical 
WVD  (fig.  8(b))  and  the  experimental  WVD  (fig.  6(d)),  show  a  large  cross  pattern  at  the  sixth  tooth  pass 
location  due  to  the  tooth  fracture.  Some  discrepancies  have  been  detected  between  the  experimental  and 
the  numerical  time  signal  at  the  fourth  and  fifth  tooth  pass  locations.  The  experimental  time  signal  con¬ 
sists  of  some  higher  frequency,  smaller  amplitude  vibration  modulation,  which  were  not  simulated.  This 
additional  modulated  signal  excited  the  14  times  rotational  speed  component,  as  shown  in  the  experimental 
frequency  spectrum  in  figure  6(d),  and,  in  turn,  is  responsible  for  the  small  differences  created  in  the  WVD. 


7 


The  FM4  parameter,  as  seen  in  figure  9(a),  shows  a  possible  reaction  as  the  pitting  started  to  oc¬ 
cur,  however,  it  does  not  provide  any  coherent  indication  of  the  severity  of  the  pitting  as  the  damage  in¬ 
creased.  In  addition,  it  does  not  provide  any  information  to  distinguish  the  pitting  on  single  or  multiple 

teeth.  ^ 

Results  from  NA4*  are  illustrated  in  figure  9(b).  It  is  obvious  from  the  figure  that  NAA  provides  a 
very  good  indication  of  the  pitting  development  on  the  pinion  tooth.  The  magnitude  of  the  parameter  in¬ 
creases  to  a  nondimensional  value  of  7  after  shutdown  no.  2  at  6.55  hr,  and  further  to  a  value  of  17  when 
the  pitting  covers  75  percent  of  the  tooth  surface  at  8.5  hr.  As  expected,  the  “locked”  denominator  in 
AA4*  provides  a  robust  indication  as  the  pitting  progresses  [Decker  1994  and  Zakrajsek  1994]. 

The  NBA*  method,  as  shown  in  figure  9(c),  displays  a  very  similar  trend  to  that  of  NAA*,  with  a 
more  robust  indication  to  the  severity  of  the  damage.  However,  both  NAA  and  NBA  did  not  provide  any 
type  of  indication  as  the  damage  spread  to  other  teeth,  and  finally  tooth  fracture. 

From  the  results,  it  is  clear  that  the  analytical  model  has  the  ability  to  accurately  model  gear  tooth 
pitting  damage.  The  WVD  method  not  only  verified  the  analytical  results,  but  also  provided  detailed  diag¬ 
nostic  information  on  the  damage.  The  WVD  method  provides  information  as  to  the  location  of  the  dam¬ 
age  on  the  gear,  and  also  the  magnitude  of  the  damage  present.  This  information,  however,  is  currently 
only  discemable  by  the  human  eye.  A  post-processing  scheme,  such  as  a  neural  network  algorithm,  is 
needed  to  process  the  data  provide  by  the  WVD  in  an  automated  fashion. 


VI.  CONCLUSIONS 

1.  A  modal  synthesis  methodology  was  developed  to  simulate  the  dynamics  of  gear  transmission 
systems.  While  the  computational  efforts  have  been  greatly  reduced  by  the  modal  transformation,  the  nu¬ 
merical  results  generated  maintain  good  accuracy. 

2.  Gear  tooth  damage  due  to  wear  and  pitting  can  be  approximated  by  amplitude  and  phase 
changes  in  the  gear  mesh  stiffness  model.  The  gear  mesh  model  can  easily  be  incorporated  into  the  global 
transmission  system  for  dynamic  simulations. 

3.  The  WVD  method  provides  a  comprehensive  representation  of  the  vibration  signal.  It  was  suc¬ 
cessfully  used  to  verify  the  analytical  model. 

4.  The  WVD  method  provides  detailed  information  that  can  be  used  to  detect  and  locate  gear 
tooth  faults. 

5.  The  NAA*  and  NBA*  parameters  show  good  reactions  to  the  initial  pitting  damage  and  very  nice 
indications  for  the  growth  and  severity  of  the  pitting  damage.  However,  their  indications  for  the  tooth 
fracture  are  unclear. 
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Figure  1  .—Schematic  of  a  rotor-gear  bearing  system. 


Figure  2. — Spiral  bevel  gear  test  rig  at  NASA  Lewis  Research  Center. 
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Figure  3. — Pictures  of  the  damaged  pinion  teeth,  (a)  5.5  hours,  (b)  6.55  hours,  (c)  8.55  hours,  (d)  1 0.03  hours. 
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Figure  4. — Pictures  of  the  damaged  pinion  teeth,  (a)  12.03  hours,  (b)  14.53  hours,  (c)  16.16  hours,  (d)  17.79  hours. 
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Figure  7. — Gear  mesh  stiffness  models  to  simulate 
damage  on  pinion  gear  teeth,  (a)  Single  tooth  damage 
model,  (b)  Multiple  tooth  damage  model. 
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Figure  8.— Time  signal,  WVD  plot,  and  frequency  spectra 
of  the  damaged  pinion  teeth  from  analytical  simulations, 
(a)  12.03  hours  (sin^  ^  tooth  damage),  (b)  17.79  hours 
(damage  on  three  consecutive  teeth). 
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